pub fn new_birthday_probability(n: u32) -> f64 {
    let mut p = 1.;

    for i in 0..n {
        p *= (365 - i) as f64 / 365 as f64;
    }

    1. - p
}
